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GENERALIZED ELECTROCHEMICAL CELL 
STATE AND PARAMETER ESTIMATOR 

FIELD OF THE INVENTION 
[0001] The present invention relates to electrochemical power 
sources, and more particularly to determining states and parameters of 
an electrochemical power source, 

5 

BACKGROUND OF THE INVENTION 
[0002] Electrical power sources, such as batteries or fuel cells, are 
made up of an array of electrochemical cells. The state of charge 
(SOC) of the electrochemical cell (EC) indicates the amount of usable 

1 0 energy stored within the EC at a given time. It is analogous to the 

amount of fuel in a fuel tank. In order to improve EC life and overall EC 
performance, the EC must be maintained at an optimal SOC. The 
state of health (SOH) of an EC indicates the physical condition of the 
EC. The SOH is based on observing particular EC parameters to 

1 5 determine the extent of unobservable, internal damage or degradation. 
[0003] An EC can be considered a system that has current and heat 
flow as inputs and responses of voltage and temperature. The voltage 
response of the EC to current stimulation is described using a 
combination of voltages and a current source. These voltages include 

20 an equilibrium potential, voltage resulting from hysteresis effects, 
voltage resulting from ohmic resistance drops, and voltage resulting 
from dynamics in the EC, or double-layer voltage. Each of these 
voltages can be described by either an algebraic function, a differential 
equation or a convolution integral. The current source is indicative of a 

25 self-discharge current of the EC. 

[0004] During use, resistance, equilibrium potential, voltage 
resulting from hysteresis effects, voltage resulting from ohmic 
resistance drops, double-layer voltage, and self-discharge current are 
not directly measurable. These values influence the SOC and the 

1 



parameters that indicate the SOH of the EC. Because they are not 
directly measurable, it is often difficult to accurately determine the SOC 
and SOH of the EC. Therefore, it is difficult to maintain the EC at an 
optimal SOC or to determine when the SOH has degraded to a point 
5 where EC performance is significantly effected. 

SUMMARY OF THE INVENTION 
[0005] Accordingly, the present invention provides a method of 
estimating a state of charge (SOC) and a state of health (SOH) of an 

10 electrochemical cell (EC). The method includes modeling the EC with 
a linear equation, measuring a terminal current of the EC, measuring a 
terminal voltage of the EC and measuring a temperature of the EC. 
The linear equation is processed through a time-varying state and 
parameter estimator based on the terminal current, the terminal voltage 

15 and the temperature to determine states and parameters of the EC. 
[0006] In one feature, the states include the SOC and ohmic 
resistance of the EC. 

[0007] In another feature, the parameters include a variable, a 
charge variable and a discharge variable of the EC. The SOH is 
20 determined based on the variable, the charge variable and the 
discharge variable. 

[0008] In another feature, the linear equation is a process model 
that models changes in the states and parameters of the EC based on 
current. 

25 [0009] In another feature, the method further includes determining 
synthesized inputs based on the terminal current, the terminal voltage 
and the temperature. Processing is further based on the synthesized 
inputs. 

[0010] In still another feature, the step of processing the linear 
30 equation includes predicting current states based on prior states and 
predicting current error of the states based on prior errors of the states. 
Current state gains are determined based on the current error of the 
states and the current states are updated based on the prior states and 
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the current state gains to provide the states of the EC. The states 
include the SOC and ohmic resistance of the EC. The state errors are 
continuously updated. 

[001 1] In yet another feature, the step of processing the linear 
5 equation includes predicting current paranneters based on prior 
parameters and predicting current error of the parameters based on 
prior errors of the parameters. Current parameter gains are 
determined based on the current error of the parameters and the 
current parameters are updated based on the prior parameters and the 

10 current parameter gains to provide the parameters of the EC. The 
parameters include a variable, a charge variable and a discharge 
variable of the EC. the parameter errors are continuously updated. 
[001 2] Further areas of applicability of the present invention will 
become apparent from the detailed description provided hereinafter. It 

15 should be understood that the detailed description and specific 

examples, while indicating the preferred embodiment of the invention, 
are intended for purposes of illustration only and are not intended to 
limit the scope of the invention. 

20 BRIEF DESCRIPTION OF THE DRAWINGS 

[0013] The present invention will become more fully understood 
from the detailed description and the accompanying drawings, wherein: 
[0014] Figure 1 is a functional block diagram of a system including a 
power source, a load powered by the battery and a generator to 

25 recharge the power source; 

[0015] Figure 2 is a schematic diagram of the power source; 
[0016] Figure 3 is a graph illustrating state of charge (SOC) versus 
open-circuit voltage (OCV) for an exemplary lead-acid (Pba) 
electrochemical cell; 

30 [0017] Figure 4 is a signal flow diagram of an electrochemical cell 
state and parameter estimator according to the present invention; 
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[0018] Figure 5 is a flowchart illustrating the electrochemical cell 
state and parameter estimation control according to the present 
invention; 

[001 9] Figure 6 is a flowchart illustrating processing of synthesis 
5 functions for the electrochemical cell state and parameter estimation 
control of Figure 4; 

[0020] Figure 7 is a flowchart illustrating state estimation for the 
electrochemical cell state and parameter estimation control of Figure 4; 
and 

1 0 [0021 ] Figure 8 is a flowchart illustrating parameter estimation for 
the electrochemical cell state and parameter estimation control of 
Figure 4. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 
15 [0022] The following description of the preferred embodiment is 
merely exemplary in nature and is in no way intended to limit the 
invention, its application, or uses. For purposes of clarity, the same 
reference numbers will be used in the drawings to identify similar 
elements. 

20 [0023] Referring now to Figure 1 , a system 10 includes a 

regenerative power source 12 that is discharged as it powers a load 
14. The power source 12 is made up of an array of electrochemical 
cells (EC*s) (not shown). The EC's are selectively recharged by a 
generator 16. Operation of the generator is controlled by a controller 

25 18. Sensors 20 communicate with the EC's and generate a terminal 
voltage signal (Vt), a terminal current signal (It) and a temperature 
signal (T). The controller 18 processes the signals according to the 
estimator control of the present invention to estimate states and 
parameters of the EC's. The controller 1 8 signals the generator 1 6 to 

30 operate based on the estimates. 

[0024] Although the system 10 is a generic power system, it is 
indicative of a vehicle power system. In such a case, the power source 
12 powers loads such as headlamps, heaters, power seats, windows, 
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and the like. The vehicle's engine drives the generator 16 (i.e. 
alternator), to produce recharging current. 

[0025] Referring now to Figure 2, the power source is schematically 
modeled in terms of an equilibrium potential or open-circuit voltage 
5 (Vo), a voltage resulting from hysteresis effects (Vhyst), a voltage 

resulting from ohmic resistance drops (Vohmic), a voltage resulting from 
dynamics in the battery or double-layer voltage (Vdi) and a self- 
discharge current (Isd). Vo can be determined based on Vdi and Vhyst. 
Vdi and Vhyst are determined using an open loop model. These 

10 voltages are subtracted from a measurable terminal voltage (Vt) to 
determine Vo- Knowing Vo is sufficient to estimate the SOC of the 
battery given a function that maps Vo to SOC. An empirical fit of the 
SOC to Vo for an exemplary power source is shown in Figure 3. The 
exemplary power source of Figure 3 is a PbA battery. Although the 

15 power source model of Figure 2 does not include all of the existing 
behaviors, testing and experience has shown that this model Is 
sufficient for SOC estimation, power control and diagnostics for 
automotive applications. 

[0026] The present invention provides a general state and 
20 parameter estimator to estimate the SOC and SOH of an EC. The 
SOC and SOH estimates are used to operate and/or diagnose the EC 
or array of EC's. The estimator is a Kalman Filter (KF) -based state 
and parameter estimator. A standard Kalman Filter (KF) combines the 
concept of recursive least squares estimation with a state space model 
25 and noise model to optimally estimate a system characteristic. The KF 
estimates the characteristic of the system based on a linear model of 
the system. The linear model of the system includes a linear dynamic 
model, a description of the covariance of the errors in the system 
dynamics, and a description of the covariance of the errors in system 
30 measurements. The linear model is expressed in state space form as: 
Xk+1 = AkXK + BkUk + Wk (Process Model); and 
yk = FkXk + GkUk + rk (Measurement Model) 

where: 
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Xk is the state vector (i.e. what is to be estimated); 
Yk is the measurement vector (known); 
Uk is the deterministic input to the battery (known); 
Wk is the process noise (i.e. random inputs into the system); 
5 rk is the measurement noise (i.e. random noise that affects 

measurements); and 
Ak, Bk, Fk, Gk are the system matrices (known). 

[0027] The KF estimates the characteristic by using a form of 
10 feedback control. In short, the KF estimates the characteristic at some 

time and then obtains feedback in the form of noisy measurements. 

The KF equations fall into two groups: time update equations 

(predictor) and measurement equations (corrector). The time update 

equations project the current characteristic and error covarlance 
15 estimates to obtain estimates for the next time step. The measurement 

update equations incorporate new measurements into the estimates to 

obtain an improved or updated estimate. 

[0028] The KF-based estimator of the present invention introduces 
a mechanism that propagates state and parameter estimates (x and p, 

20 respectively) based on EC dynamics. The estimator propagates the 
estimate error covariance based on EC dynamics and knowledge of 
how noise affects the EC. The EC is a non-linear system. However, 
the KF-based state and parameter estimator simplifies the system by 
either linearizing non-linearities or removing them. Further, the KF- 

25 based state and parameter estimator accounts for time varying 

parameters to more quickly improve estimate errors. The controller 
processes the KF-based state and parameter estimator and controls 
other components of the system 10 in response to the state and 
parameter estimates. 

30 [0029] The state estimates generally include SOC, ohmic resistance 
(Ro) and double layer voltage (Vdi). As used herein, the acronym VRV 
indicates the state estimates. The parameter estimates include a 
variable (f), a charge variable (bcharge) and a discharge variable 
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10 



15 



(bdischarge)* As used herein, the acronym fbb indicates the parameter 
estimates. The variables f , bcharge and bdischarge are based on the 
dynamic impedance (Zdyn) of the EC. Zdyn is based on the resistance 
(Ret), double layer capacitance (Cdi) and time constant (t), 
[0030] As used herein, the notation k indicates the present Iteration 
of time. The notation k|k indicates the present iteration evaluated using 
information presently available (i.e. updating or correcting). The 
notation k(k-1 indicates the present iteration evaluated using past 
(stored) information (i.e. estimating). Additionally, the notation z""^, 
more specifically 2"\ indicates a unit delay. 
[0031] In general, the ECSP estimator includes six processing 
models: a synthesis model, a state prediction model, a parameter 
prediction model, a measurement model, an estimate correction model 
and an error propagation model. The notation {^) indicates variables 
determined using the measurement model. Initially, Vt, It and T are 
measured for the EC. These values are processed through the 
synthesis model to provide an input vector Usyn- The input vector is 
provided as: 

h 



^ choree, k 



disch2Sge,k 



20 



where: 



^ charge = max(0, /) 

^rf«cAarg* =niin(0,/) 



25 [0032] The measurement model is provided as: 
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^ ' ~ fqtov (Qsoc,Jt|it-i » ^)"^ 1000 
-1 



R 



o,k\k-\ ~^ ^dLk\k-\ 



''r.*|*-l 



1000, 



\+b 



1000 



y = 



A 



[0033] The state prediction model propagates (predicts) the 
estimated state and parameter values (Xk,Pk) and the covariance (P) of 
the respective estimate errors based on prior information. The state 
prediction model is provided as: 

VCE " fchargeff {QsOC,k\k 

QsocMf^-i ~ Qsoc.k\k {^CE * ^rf ) * 2500 

^dLklk-l ~ fk\k '^dl,k\k '^K,k\k 1000 ^ 



1000 



QsOC,k\k-\ 
^dKk\k-l 



The parameter prediction model is provided as: 







' ^o.k\k 


" 


" fk\k 






Kk]fc 


^d,k\k-\ 




^d,k\k 

Ro,k\k-l 


Pk\k-\ ' 




fk\k-\ 
^c,k\k-l 



[0034] The next measurement is obtained and the estimator 
updates (corrects) the predicted Xk and pk based on the current 
measurements. The predicted covariance (P) of the respective 
estimate errors is updated based on the measurements. The order of 
these operations can vary. 

[0035] The estimator implements the following relationships and 
matrices for state estimation: 
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c = 



1 0 
0 1 
0 0 



0 
0 



0 a:^,.,„.,/1000 
0 0 

fviQ (GsOC.*|* ' ) ^^^^ 1 



1 

0 



^VRV 



[0036] The state estimate correction model includes a state gain 
matrix (Kvrv) based on the state error covariance (PvRv.k|k-i), the matrix 
CvRv and Rvrv- Kvrv is calculated according to the following equation: 

^VRV ~ ^VRV,k\k-\ ' ^VRV ' {^VRV ' ^VRV ,k\k-\ ' ^VRV ^VRV ) 



[0037] A previously stored state estimate (Xk|k-i) is updated to 
determine a present state estimate (Xk|k). The current state estimate is 
determined based on Xk|k-i, Kvrv and an a posteriori error (evRv). The 
10 error is defined as the difference between the actual Vt and the terminal 

A 

voltage based on the measurement {Vt ). The present state estimate is 
determined according to the following equation: 

^VRV,k\k " ^VRV,k\k-l ^VRV ' ^VRV 



15 [0038] The state estimate is propagated based on the present state 
estimate according to the following equation: 

^VRVMi\k ~ ^VRV ' ^VRV,k\k ^VRV ' ^VRV 

[0039] The state estimate error covariance (Pvrv) is propagated 
based on the updated state estimate error covariance PvRv.k according 
20 to the following equation: 

^VRVMl\k ~ ^VRV ' ^VRV,k\k ' ^VRV QvRV 
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[0040] To summarize state estimation, the estimator inputs include: 
Vt,k. It.k, Tk, Pfbb.k|k-i, xvRv.k|k-i and PvRv.k|k-i. The values Vt,k, lt,k and Tk 
are directly measured from the EC. The values xvRv,k|k-i and PvRv,k|k-i 
are previously determined and retrieved from memory. The estimator 
5 outputs include: xvRv,k+i|k, PvRv.k+i|k, xvRv.k|k and PvRv.k|k. The values 
XvRv.k+i|k and PvRv.k+i|k are stored in memory for use during the next 
iteration. The value xvRv.Kjk defines the estimated SOC, ohmic 
resistance (Ro) and double layer voltage (Vdi). These values are used 
to operate the system and/or diagnose the EC. 
10 [0041] The estimator implements the following relationships and 
matrices for parameter estimation: 

^rf/.Jfc-l|Jfc-l ~ ^VRV^k-\\k-l 



fk-\ ~ Pjbb,k-\\k- 



^Jbb ~ 



dl,k-l\k-l ^QQQ JpQQ 

Qjbb ~ f/bb-Q^'"^ 

^Jbb = (l + fk-\ ) ' PvRV,k\k 



[0042] To estimate pfbb, the estimate correction model includes a 
15 parameter gain matrix (Kfbb) based on the parameter error covariance 
(Pfbb,k|k-i), the matrix Cfbb and Rfbb. K^b is calculated according to the 
following equation: 

^fbb ~ ^Jbb,k\k-l ' ^fbb ' {^Jbb ' Pfbbk\k-l ' ^fbb ^Jbb ) 

20 [0043] A previously stored parameter estimate (Pk|k-i) is updated to 
determine a current parameter estimate (pk|k). The current parameter 
estimate is determined based on pk|k-i, Kfbb and an a posteriori error 
(efbb,k). The error is defined as the difference between the actual Vdi 
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and the Vm , which is based on the measurement. The current 
parameter estimate is determined according to the following equation: 

Pjbb,k\k ~ Pjbb.kVc-i ^Jbb '^Jbb 

[0044] The parameter estimate error covariance (Pfbb) is propagated 
according to the following equation: 

PjbbM\\lt ~ ^Jbb ' ^Jbb.k\k ' ^Jbb + Qjbb 

[0045] The parameter estimate Is propagated based on the present 
state estimate according to the following equation: 

PjbbM\\fi ~ ^fbb ' PjbbMk ^Jbb ' "jbb 

[0046] To summarize parameter estimation, the estimator inputs 
include: lt,k, Tk, Pfbb.k|k-i, xypv.kik, xvRv.k-i|k-i, PvRv.k|k, Xfbb.k|k-i and Pfbb.k|k-i. 
The values lt,k and Tk are directly measured from the EC. The values 
pfbb.k|k-i, xvRv,k|k-i and PvRv.k|k-i are previously detemiined and retrieved 
from memory. The estimator outputs include: Pfbb,k+i|k, Pfbb,k+i|k, Pmb.kik 
and Pfbb,k|k- The values Pfbb.k+i|k and Pfbb,k+i|k are stored in memory for 
use during the next iteration. The value Pfbb,k|k defines the estimated 
ohmic resistance (Ro), f, bcharge and bdischarge- These values are used to 
operate the system and/or diagnose the EC. 
[0047] Referring now to Figure 4, a signal flow for the general EC 
state and parameter (ECSP) estimator will be described. Voltage (V), 
current (I) and temperature (T) are measured for the EC at 100. At 
102, the V, I and T signals are processed through the synthesis 
functions by the controller 18. As described above, the synthesis 
functions produce the input vector signal (Usyn). The input vector signal 
(Usyn) is processed in the various models and functions of the ECSP 
estimator. 

[0048] The input vector signal is processed through the state 
predictor model to produce a state estimate signal (Xk|k-i) at 104. The 
state estimate signal is stored in a buffer at 106 for later processing. At 
1 08, the parameter predictor model processes the input vector signal to 
produce a parameter estimate signal (Pk|k-i). The parameter estimate 
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signal is stored in a buffer at 1 10 for later processing. The input vector 
signal is processed through the error predictor model at 1 12 to produce 
an error estimate signal (Pk|k-i). The error estimate signal is stored in a 
buffer for later processing at 1 14. 

5 [0049] The state and parameter estimate signals (Xk|k-i, Pk|k-i) are 
pushed from their respective buffers to the measurement model at 1 16. 
The input signal (Usyn) and corrected state, parameter and error signals 
(Xk|k, Pk|k, Pk|k, respectively) are processed in the measurement model 
with the estimated state and parameter estimate signals (Xk|k-i, Pk|k-i). 

10 The measurement model produces a voltage signal (y). 

[0050] At 1 18, the voltage signal (y) is processed through the 
correction function with the input signal (Usyn), the state and parameter 
estimate signals (Xk|k.i, Pk|k-i) and the error estimate signal (Pkik i). The 
correction function produces the corrected state signal (Xk|k), the 

15 corrected parameter signal (pk|k), and a corrected error signal (Pk|k). 
The corrected state signal is processed with the input signal through 
the state predictor model to produce the state estimate signal. The 
corrected parameter signal is processed with the input signal through 
the parameter predictor model to produce the parameter estimate 

20 signal. The corrected error signal is processed with the input signal 
through the error predictor model to produce the error estimate signal. 
[0051] Referring now to Figure 5, the ECSP estimator control will be 
described in detail. In step 500, control determines whether the system 
is on. If so, control continues in step 502. If not, control loops back. 

25 Control determines initial values for Xvrv, Ptbb, Pvrv and Pfbb in step 

502. Initial values for the power source states and parameters are pre- 
programmed in memory. Certain assumptions are required to 
determine these initial values. The initial values are selected to 
minimize the error between the actual states and the assumed states 

30 to enable quicker error reduction. Assumptions are also required for 
the measurement and process noises. In step 504, control sets a 
counter (k) equal to 1 . 
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[0052] In step 506, the voltage(Vt), current (h) and temperature (T) 
are measured. In step 508, control processes the measurements with 
the synthesis functions to produce the input signal. This is described in 
further detail with reference to Figure 6. Control determines the 
5 corrected state estimates and the corrected parameter estimates in 
steps 510 and 512, respectively. Calculation of the corrected state 
estimates and the corrected parameter estimates is described in more 
detail with reference to Figures 7 and 8, respectively. 
[0053] In step 51 4, control operates and/or diagnoses the system 

10 10 based on the corrected state and parameter estimates. In step 516, 
control determines whether the system is off. If not, control increments 
the counter in step 518 and loops back to step 506. If so, control ends. 
[0054] Referring now to Figure 6, control determines the input signal 
(Usyn) based on the measurements Vt, It and T in step 600. In 

15 particular, control determines I, Ic and Id. As described above, the input 
signal (Usyn) is a vector including Vt, I, Ic, Id and T, 
[0055] Referring now to Figure 7, the process to determine the state 
estimate will be described. In step 700, control loads the matrices 
used to determine the VRV state estimate. These matrices are shown 

20 above. Control calculates the state gain matrix (Kvrv) in step 702. In 
step 704, control estimates the current state estimate (xvRv.k|k-i) based 
on prior information. In step 706, control estimates the current error 
covariance estimate (PvRv.k|k-i) based on prior information. In step 706, 
control updates the stored state estimate (xvRv.k|k). In step 710, control 

25 updates the stored error covariance estimate (PvRv.k|k). In step 71 2, 
control stores the updated estimates XvRv,k|k and PvRv.kik- 
[0056] Referring now to Figure 8, the process to determine the 
parameter estimate will be described. In step 800, control loads the 
matrices used to determine the parameter estimate (fbb). These 

30 matrices include are described above. Control calculates the 

parameter gain matrix (Kfbb) in step 802. In step 804, control estimates 
the current parameter estimate (Pfbb.k|k-i) based on prior information. In 
step 806, control estimates the current error covariance estimate 
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10 



15 



20 



25 



(Pfbb,k|k-i) based on prior information, in step 810, control updates tlie 
stored parameter estimate (ptbb.kik)- In step 810, control updates the 
stored error covariance estimate (Pfbb.kik). In step 812, control stores 
the updated estimates p^b.kik and Pfbb.kik. 

[0057] In sum, the ECSP estimator process Vt.k, lt,k and T through 
the various models to provide estimate stats and parameter outputs. 
The state outputs are provided as: 



Qsoc 

R 



The parameter outputs are provided as: 



PfbbMl\k " 



f 



^oMisch dirge 



The state and parameter outputs are used to operate and diagnose the 
power source. 

[0058] The present invention further provides variants of the general 
ECSP estimator described in detail above. The estimator variants 
described herein operate using the same logic as described in Figures 
5 through 8. However, the relationships and matrices that power the 
variants vary based on the particular design goal of each. 
[0059] A first variant is directed toward symmetric or asymmetric 
impedance batteries. The first variant uses a general characterization 
of the EC, introduces hysteresis effects and uses multiple independent 
descriptions of estimate errors. The first variant estimator can be used 
on either symmetric or asymmetric impedance EC's with only 
calibration changes. The use of characterization data to improve 
estimator performance with little additional processing. The 
characterization data includes the rate of change on R, f and b with 
respect to state variables such as Q and T. Exemplary 
characterization data includes: dR/DQ, dR/dT, df/dQ, df/dT, db/dQ and 
db/dT. 
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[0060] The first variant inputs are provided as: 



ch arg e 



dischsarge _ 



' syn 



= [vj 



where the synthesis model provides: 



^ ~" ^SD [Qsoc,k\k »^) 

= max(/,0) 
I, =min(/,0) 



1 0 [0061] The equations of the state prediction model for the first 
estimator variant are provided as: 

VcE " fchargeeff [QsOC,k\k 

fo ~ ff (Csociti* 
= f^QsocMk^T) 

QsocM^i-i ~ Qsoc,k\k (^c£ ' K ) * 2600 

^rf/.*|A-l ~ {f offset. k\k fo)'^dl,k\k ^ offset, k\i(^ ) ' ^JQQQ ^ ' \QQQ j 



Q 



SOC,k\k-\ 



dl,k\k-\ 
QH,k\k-l 



[0062] The equations of the parameter prediction model for the first 
1 5 estimator variant are provided as: 
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dRdQ = f,^{Qsoc.k\k'T) 

SOC,ik|t»' / 

dfdQ = /rf^rf, (Gsoc.*!* 

/<!ffiter,*|*-l ~ foffsei.k\k ''S^G * 

dbdQ = f^^a,[Qsoc.k\k^T) 
dbdT = f^^(Qsoc,kik^T) 



QsOC.k\k ' ^ ) ^ 

Ts 



+ dRdT- 



+ dfdT- 



Kff.et.k\k-i = + dbdQ 



Ts 



-^-dbdT 



Ts 



( t-(1-z-') ^ 
Ts 



f offset, k\k-l 
offset. k\k-\ 



[0063] The equations of the measurement model for the first 
estimator variant are provided as: 

^/ /q/oV (2sOC,A|Jfc-l » ^) + JQQQ ' ^ offset, k\k-\ ) ^rf/,it|*-l foutoVf, ^SOC,k\k-\ ' 2/y 

K// = (/o + fofset,k\k'-\)' ^dUk\k O"*" (^o ^qflr5e/,ife|it-l )*^2QQQ * ^ JQQQ 



[0064] The equations of the estimate correction model for the first 
estimator variant are provided as: 
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^VRV \fdVdQ (QsOC .Jk|Jk-l > ^ ) fdRdQ (QsOC,Jt|Jt-l > ^ ) * ^ ^ 1 fdVdQf, (CsOOtl* -P ^ )J 



SOC,it|jt-l 



offset, k\k'-\ 



-VTiV 



SOC^\k 



qffse:,k\k 



VRV ^^VRV 



" ^VRV,k\k-\ 

^VRV,k\k ~ ~ ^V/?V * ^VRV )' ^V7?V.jt|ft-l 
'^ = /r(GsOC.it!JtT) 

~ (l + fplus)' ^VRV,k\k 



^fbb^klk-l 



f offset, k\k-l 
offset, k\k-l _ 



^y&i? "~ ^dl,k\k 



^fbb ^fbb,k\k-~\ * * (C'ybfr • ^fbb.k\k-\ ' ^fbb ^fbb ) 



^fbbMk 



b 



offset, k\k 



^fbb,k\k-l ^y&ft ' ^fbb 



^fbb,k\k ~ ~ ^ybt ' ^fbb )* ^yb/?,JtIit-l 

[0065] The equations of the error propagation model for the first 
estimator variant are provided as: 
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10 



[0066] 



fo ~" // {QsOC.k\k ' 

"1 0 0 0' 

0 1 0 0 

0 0 + 0 

0 0 0 1 

QvRV ~ /qvrv {QsOC.k\/t ^) 
Qjbf, ~ /gja^ {QsOC,k]k »^» ^) 

'vRV,it)Jt-l ~ ^VRV ' ^VRV,k\k ' A//?V QvRV 
^fbb,k\k-\ ~ ^jbb,k\^ '^Qjbb 
^k\k'i ^ {^VRV,Jtjit-l ' ^Jbb,k\k-l } 

The outputs of the second variant are provided as: 



^k\k " 



Pk\k 



Q 



SOC 

^offset 
f offset 
^offset 



^k\k l^VRV ' ^fbb } 



[0067] A second variant learns the Impulse response of the 
linearized system to learn the impedance of the system. The second 
variant uses a finite impulse response and ordinary differential equation 
parameters to adapt to the power source characteristics using a single, 
unified description of the estimate errors (i.e. single matrix). The 
impulse response is learned as a set of discrete values that are spaced 
in time. The difference between the predicted output and the 
measured output is used to drive updates to the learned impulse 
response. 

[0068] The input vector of the second variant is provided as: 



^syn = 



15 



18 



where: 

^ ~ ^1 ~ ^ SD{Qsoc.k\k »^) 

Icins, =niax(/,0) 
hjns, =min(/,0) 

^ c [^r.m.« ^ cinst ^ ^ cinst ^ 

^ d ~~ djnst ^ djnst ' ^ ^djnst 



[0069] The equations of the state prediction model for the second 
estimator variant are provided as: 

VcE ~ fchvurgeeff {QsOC,k\k 

r = fr(QsocM^^^) 

fo ~ ff {QsOC,k\f( ' ^) 

= fb(Qsoc,k\k^'^) 



Q 



SOC.*(*-I 



3600 



T)-sign{l)- 



^dl,k\k~l 



^k\k-l " 



^dl,k\k'l 



1000 



1000 



[0070] The equations of the parameter prediction model for the 
10 second estimator variant are provided as: 





^ck\k 


1 " 


= 7 

^d,k\k 




^Cklk-l 




_^d,k\k'\^ 



[0071] The equations of the measurement model for the second 
estimator variant are provided as: 
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^/ fatoV {.QsOC^k-l '^)"'" Qh K ' ^c,k\k-l ^rf * ^d,k\k-l K/j/i 



[0072] The equations of the estimate correction model for the 
second estimator variant are provided as: 

^ ~ IfdVdQ {QsOC.k\k'l ' ^) fdVdQii (GsOC,Jt|Jk-l » ^) ^ ^rf ] 



R = 


















-I ~ 
















e = 


V,- 


A 

V, 



^k\k ~~ 



^ offset, k\k 



' d\,k\k 
QH,k\k 

Pk\k =(l"^*^)"^fc|Jt-l 

5 [0073] The equations of the error propagation model for the second 
estimator variant are provided as: 

"1 0 0 0 0^ 

0 0 0 0 

A= 0 0 1 0 0 

0 0 0 1 0 

0 0 0 0 1^ 

[0074] The output of the second variant is provided as: 

Pk\k = 



■'c,k\k 



^d,k\k 
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[0075] A third variant is directed toward simple impedance 
estimation. In cases such as diagnostics (e.g. diagnosing weak EC's) 
or determining the characteristics of chemical reactions occurring 
within an EC (e.g. interface characteristics of a fuel cell), it is only 
necessary to learn the resistance and voltage dynamics of the EC. 
The third variant achieves this using a known open-circuit voltage (Vo). 
[0076] The input vector of the third variant is provided as: 



^syn = 



10 where: 



= max(/,0) 
I, =min(/,0) 



[0077] The equations of the state prediction model for the third 
1 5 estimator variant are provided as: 

^k\k-l ~ IVdhklk-l ] 



[0078] The equations of the parameter prediction model for the third 
estimator variant are provided as: 
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fk\k-l ~ fk\k 



Pk\k-\ " 



R 



fk\k-\ 
^k\k-\ 



[0079] The equations of the measurement model for the third 
estimator variant are provided as: 



'•=/r(«2so«|it'^) 

^rf/ = fk\k-\ ' ^dLk\k ' )+ ^k\k-l 



1000 



1000 



3^ = 



5 [0080] The equations of the estimate correction model for the third 
estimator variant are provided as: 
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=[/ 1] 



R 



o,k\k-l 



^VRV 



^VRV ~ ^ 



VRV,k\k-\ 



^VRV ' ipVRV ' ^\ 



^VRV,k\k 



dl.k\k 



VRV^k\k 



~ ^VRV^klk-X ~^ ^VRV '^VRV 



-1 * ^VRV ^VRV ) 



~ fk[k )' fvRV,k\k 

fk\k-\ 



^fbb,k\k~l 



^fbb ^dl,k\k 



'k\k-\ _ 

A 



^JbbM^ 



^Sbb,k\k ~ ^Jbb * )* ''^ybi^,*!/:-! 

[0081] The equations of the error propagation model for the third 
estimator variant are provided as: 

^ri 0 ' 

Krv - Q f 

QvRV ~ /fiv^Kv i^SOCMf^ ? ^» ^) 
Qfbb ^ ^SOCM\k '^»^) 



VRW ,k\k- 



~ ^VRV ' ^WRV ,k\k * ^ QwRV 



^ y&t,Jk|Jt-l ' Sbb,k^ Qjbb 



[0082] The output of the third variant is provided as: 
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Pk\k = 



R 



o,k\k 
fk\k 

bk\k . 



^k\k " '/v/?V,t|jt ' ^Jbb,k\k } 

[0083] EC's are significantly affected by temperature and interact 
with any heat exchange system by generating and absorbing heat 
during reactions. A fourth variant provides a way to include 
thermodynamics and heat transfer to estimate the temperature of the 
reaction. 

[0084] The input vector of the fourth variant is provided as: 

/ 



^syn = 



celLmeas 

T 

„ coolant 



10 where: 



= max(/,0) 
=min(/,0) 



[0085] The equations of the state prediction model for the fourth 
1 5 estimator variant are provided as: 
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VCE " fchM^eeff ^SOC.k\k '^re//.A:|* ) 



3600 



1000 



1000 



mC 



p^cell 



Xk]}c-\ = 



SOC,k\^-\ 
^dl,k\k-\ 



' re//,*|i5:-l 



[0086] The equations of the parameter prediction model for the 
fourth estimator variant are provided as: 

fk\k-\ ~ fk\k 



Pk\k-\ ~ 



fk\k-l 

K,k\k-\ 



[0087] The equations of the measurement model for the fourth 
estimator variant are provided as: 



- fQtov{Qsoc.k\k-i^'^)~^ 1000 



k\k 



+ b. 



1000 



' \+b. 



d,k\k- 



Uooo 



Tcell — T^eU,k\t-l 



y = 



V, 

A 
A 

Tcell 



[0088] The equations of the estimate correction model for the fourth 
estimator variant are provided as: 
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„ _ /rfVrfe(2soc.*i*-i'^) 

0 0 0 1 



CsOC,k\k- 



R 



Pfbb,k\k-\ 



o,*|Jt-l 
^dl,k\k-l 
_ '^cell,k\k-\ 



fk\k-l 
^c,k\k-\ 
^d,k\k-\ _ 



WRV,k\k 









T 

cell, meets _ 




^VRV,k]}c-'\ 










^o,k\k 






^dl,k\k 







V, 

A 

T cell 



VRV ' ipVRV ' ^VRV,k\k-l ' ^VRV ^VRV ) 



•^V«V,ik|*-l ^VRV ' ^VRV 



^VRV,k\k *" ^V/?V ■ C'vRV )' ^\ 



VRV,k\k 

A 

^y&fr ^ ^dl,k\k dl 



VRV,k\k-\ 



^fbb = 


Pjbb,k\k-l 


^%b 'V-'/bb ' ^fbb.klk-'l ' ^fbb ^Jbb 






fk\k 




Pjbb,k\k 




^c,k\k 


— Pfbb,k\k-l ^Jbb ' ^Jbb 






J^d,k\k _ 




^fbb,k\k 


= \ 


/ ^/bb ' ^fbb)' ^fbb*k\k-\ 



[0089] The equations of the error propagation model for the 
estimator variant are provided as: 
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1 


0 


0 


0 




0 


1 


0 


0 


A'RV ~ 


0 


0 


/tl* 


0 




0 


0 


0 















QvRV /gv»v iQsOC,k\k »^»^) 
^VRV,k\k~l ~ ^V/?V ■ ^VRV,k\k ' ^VRV 2l 



VRV 



Jbb,k\k- 



[0090] The output of the fourth variant is provided as: 



CSOC.k\k 
^dUk\k 



Pk\k = 



^cell 



R 



fk\k 



^c.k\k 
J^d,k\k _ 



p 



} 



[0091] It is also anticipated that an array of ECSP estimators can be 
5 used to estimate multiple EC's in an array of EC's (e.g. battery). In 
another instance, a single estimator can operate using the average, 
minimum or maximum measurements of an array of EC's. In still 
another instance, a combination of estimators is used either chained 
together or working on different averages. For example, one estimator 

10 can operate on the average voltage of an array of EC's and estimate 
the average open circuit voltage (Vq) of the array. The average Vo is 
then used by a bank of estimators (e.g. third variant type) that 
determine impedance for the individual EC's. In yet another instance, 
a single estimator is fed by multiple estimators. In a fuel cell, for 

15 example, multiple estimators estimate the Vo of a individual EC's 
making up the fuel cell based on temperature. The Vo values are 
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processed by an impedance estimator to determine properties of the 
chemical reactions within the fuel cell. 
[0092] Those skilled in the art can now appreciate from the 
foregoing description that the broad teachings of the present invention 
5 can be implemented in a variety of forms. Therefore, while this 

invention has been described in connection with particular examples 
thereof, the true scope of the invention should not be so limited since 
other modifications will become apparent to the skilled practitioner 
upon a study of the drawings, the specification and the following 
10 claims. 
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